Kelpie: generating full-length 'amplicons' from whole-metagenome datasets.

INTRODUCTION: Whole-metagenome sequencing can be a rich source of information about the structure and function of entire metagenomic communities, but getting accurate and reliable results from these datasets can be challenging. Analysis of these datasets is founded on the mapping of sequencing reads onto known genomic regions from known organisms, but short reads will often map equally well to multiple regions, and to multiple reference organisms. Assembling metagenomic datasets prior to mapping can generate much longer and more precisely mappable sequences but the presence of closely related organisms and highly conserved regions makes metagenomic assembly challenging, and some regions of particular interest can assemble poorly. One solution to these problems is to use specialised tools, such as Kelpie, that can accurately extract and assemble full-length sequences for defined genomic regions from whole-metagenome datasets. METHODS: Kelpie is a kMer-based tool that generates full-length amplicon-like sequences from whole-metagenome datasets. It takes a pair of primer sequences and a set of metagenomic reads, and uses a combination of kMer filtering, error correction and assembly techniques to construct sets of full-length inter-primer sequences. RESULTS: The effectiveness of Kelpie is demonstrated here through the extraction and assembly of full-length ribosomal marker gene regions, as this allows comparisons with conventional amplicon sequencing and published metagenomic benchmarks. The results show that the Kelpie-generated sequences and community profiles closely match those produced by amplicon sequencing, down to low abundance levels, and running Kelpie on the synthetic CAMI metagenomic benchmarking datasets shows similar high levels of both precision and recall. CONCLUSIONS: Kelpie can be thought of as being somewhat like an in-silico PCR tool, taking a primer pair and producing the resulting 'amplicons' from a whole-metagenome dataset. Marker regions from the 16S rRNA gene were used here as an example because this allowed the overall accuracy of Kelpie to be evaluated through comparisons with other datasets, approaches and benchmarks. Kelpie is not limited to this application though, and can be used to extract and assemble any genomic region present in a whole metagenome dataset, as long as it is bound by a pairs of highly conserved primer sequences.

Gamarada debralockiae gen. nov. sp. nov.-the genome of the most widespread Australian ericoid mycorrhizal fungus.

This study describes a novel ericoid mycorrhizal fungus (ErMF), Gamarada debralockiae Midgley and Tran-Dinh gen. nov. sp. nov. Additionally, catabolism was explored from a genomic perspective. The nuclear and mitochondrial genomes of G. debralockiae were sequenced. Morphological characteristics were assessed on various media. Catabolic genes of G. debralockiae were explored using SignalP and dbCAN. Phylogenetic comparisons were undertaken using Phylogeny.fr. The 58.5-Mbp draft genome of G. debralockiae contained 17,075 putative genes. The complete mitochondrial genome was 28,168 bp in length. In culture, G. debralockiae produces slow-growing non-sporulating colonies. Gamarada debralockiae has many putative secreted catabolic enzymes. Phylogeny indicated G. debralockiae was distinct from known ascomycetous ErMF: Pezoloma ericae, Meliniomyces spp., Oidiodendron spp., and Cairneyella variabilis. It is closely related to many undescribed plant root-associated fungi and its nearest described relative is Hyphodiscus brevicollaris. Gamarada debralockiae has been recovered from virtually all Australian ericoid mycorrhizal studies and biogeographic data suggests the taxon is widespread in Australia. Gamarada debralockiae has similar catabolic potential to C. variabilis and co-occurs with C. variabilis at Australian sites. Plants that host multiple ErMF may benefit from subtle differences in catabolism that improve access to nitrogen and phosphorus from within recalcitrant organic matter.

First evidence of Pezoloma ericae in Australia: using the Biomes of Australia Soil Environments (BASE) to explore the Australian phylogeography of known ericoid mycorrhizal and root-associated fungi.

The prominent ericoid mycorrhizal fungus, Pezoloma ericae, has not been found in Australia to date. In the present study, internal transcribed spacer (ITS) data from the Biomes of Australia Soil Environments (BASE) was searched for evidence of P. ericae and other known ericoid mycorrhizal and root-associated taxa. ITS sequences with high identity to P. ericae, Meliniomyces bicolor, Meliniomyces variabilis, Cairneyella sp. 2, Cadophora finlandica and Woollsia mycorrhizal fungus VI were identified and their distribution in Australia visualised. This is the first evidence that P. ericae, M. bicolor and M. variabilis very likely occur on the Australian continent and provides a set of locations from which to seek isolates for further characterisation. The presence of P. ericae in South America, South Africa, and now Australia suggests a broad and ancient Gondwanan distribution for this well-studied species.

Genomic insights into the carbohydrate catabolism of Cairneyella variabilis gen. nov. sp. nov., the first reports from a genome of an ericoid mycorrhizal fungus from the southern hemisphere.

This paper describes a novel species of ericoid mycorrhizal fungus from Australia, Cairneyella variabilis, Midgley and Tran-Dinh, gen. nov. sp. nov. The genome of C. variabilis was sequenced and a draft genome assembled. The draft genome of C. variabilis is 52.4 Mbp in length, and to our knowledge, this is the first study to present a genome of an ericoid mycorrhizal fungus from the southern hemisphere. Using the SignalP and dbCAN bioinformatic pipelines, a study of the catabolic potential of C. variabilis was undertaken and showed genes for an array of degradative enzymes, most of which appear to be secreted from the hyphae, to access a suite of different carbon sources. Isolates of C. variabilis have been previously shown to utilise cellulose, carboxymethyl cellulose (CMC), cellobiose, xylan, pectin, starch and tannic acid for growth, and in the current study, putative enzymes for these processes were revealed. These enzymes likely play key roles in nutrient cycling and other edaphic processes in heathland environments. ITS phylogenetic analyses showed C. variabilis to be distinct from the fungi of the "Hymenoscyphus ericae aggregate".

Fungal identification using a Bayesian classifier and the Warcup training set of internal transcribed spacer sequences.

Fungi are key organisms in many ecological processes and communities. Rapid and low cost surveys of the fungal members of a community can be undertaken by isolating and sequencing a taxonomically informative genomic region, such as the ITS (internal transcribed spacer), from DNA extracted from a metagenomic sample, and then classifying these sequences to determine which organisms are present. This paper announces the availability of the Warcup ITS training set and shows how it can be used with the Ribosomal Database Project (RDP) Bayesian Classifier to rapidly and accurately identify fungi using ITS sequences. The classifications can be down to species level and use conventional literature-based mycological nomenclature and taxonomic assignments.

Draft Genome Sequence of Clostridium beijerinckii Ne1, Clostridia from an Enrichment Culture Obtained from Australian Subterranean Termite, Nasutitermes exitiosus.

The draft genome of Clostridium beijerinckii strain Ne1 was reconstructed from the metagenomic sequence of a mixed-microbial consortium that produced commercially significant quantities of hydrogen from xylan as a sole feedstock. The organism possesses relatively limited hemicellulolytic capacity and likely requires the action of other organisms to completely degrade xylan.

Draft Genome Sequence of Ruminoclostridium sp. Ne3, Clostridia from an Enrichment Culture Obtained from Australian Subterranean Termite, Nasutitermes exitiosus.

The draft genome sequence of Ruminoclostridium sp. Ne3 was reconstructed from the metagenome of a hydrogenogenic microbial consortium growing on xylan. The organism is likely the primary hemicellulose degrader within the consortium.

Draft Genome Sequence of Clostridium sp. Ne2, Clostridia from an Enrichment Culture Obtained from Australian Subterranean Termite, Nasutitermes exitiosus.

The draft genome sequence of Clostridium sp. Ne2 was reconstructed from a metagenome of a hydrogenogenic microbial consortium. The organism is most closely related to Clostridium magnum and is a strict anaerobe that is predicted to ferment a range of simple sugars.

Genomic and chemical insights into biosurfactant production by the mangrove-derived strain Bacillus safensis CCMA-560.

Many Bacillus species can produce biosurfactant, although most of the studies on lipopeptide production by this genus have been focused on Bacillus subtilis. Surfactants are broadly used in pharmaceutical, food and petroleum industry, and biological surfactant shows some advantages over the chemical surfactants, such as less toxicity, production from renewable, cheaper feedstocks and development of novel recombinant hyperproducer strains. This study is aimed to unveil the biosurfactant metabolic pathway and chemical composition in Bacillus safensis strain CCMA-560. The whole genome of the CCMA-560 strain was previously sequenced, and with the aid of bioinformatics tools, its biosurfactant metabolic pathway was compared to other pathways of closely related species. Fourier transform infrared (FTIR) and high-resolution TOF mass spectrometry (MS) were used to characterize the biosurfactant molecule. B. safensis CCMA-560 metabolic pathway is similar to other Bacillus species; however, some differences in amino acid incorporation were observed, and chemical analyses corroborated the genetic results. The strain CCMA-560 harbours two genes flanked by srfAC and srfAD not present in other Bacillus spp., which can be involved in the production of the analogue gramicidin. FTIR and MS showed that B. safensis CCMA-560 produces a mixture of at least four lipopeptides with seven amino acids incorporated and a fatty acid chain with 14 carbons, which makes this molecule similar to the biosurfactant of Bacillus pumilus, namely, pumilacidin. This is the first report on the biosurfactant production by B. safensis, encompassing the investigation of the metabolic pathway and chemical characterization of the biosurfactant molecule.

Multilocus variable-number tandem-repeat analysis of clinical isolates of Aspergillus flavus from Iran reveals the first cases of Aspergillus minisclerotigenes associated with human infection.

BACKGROUND: Aspergillus flavus is intensively studied for its role in infecting crop plants and contaminating produce with aflatoxin, but its role as a human pathogen is less well understood. In parts of the Middle East and India, A. flavus surpasses A. fumigatus as a cause of invasive aspergillosis and is a significant cause of cutaneous, sinus, nasal and nail infections. METHODS: A collection of 45 clinical and 10 environmental A. flavus isolates from Iran were analysed using Variable-Number Tandem-Repeat (VNTR) markers with MICROSAT and goeBURST to determine their genetic diversity and their relatedness to clinical and environmental A. flavus isolates from Australia. Phylogeny was assessed using partial ß-tubulin and calmodulin gene sequencing, and mating type was determined by PCR. Antifungal susceptibility testing was performed on selected isolates using a reference microbroth dilution method. RESULTS: There was considerable diversity in the A. flavus collection, with no segregation on goeBURST networks according to source or geographic location. Three Iranian isolates, two from sinus infections and one from a paranasal infection grouped with Aspergillus minisclerotigenes, and all produced B and G aflatoxin. Phylogenic analysis using partial ß-tubulin and calmodulin sequencing confirmed two of these as A. minisclerotigenes, while the third could not be differentiated from A. flavus and related species within Aspergillus section flavi. Based on epidemiological cut-off values, the A. minisclerotigens and A. flavus isolates tested were susceptible to commonly used antifungal drugs. CONCLUSIONS: This is the first report of human infection due to A. minisclerotigenes, and it raises the possiblity that other species within Aspergillus section flavi may also cause clinical disease. Clinical isolates of A. flavus from Iran are not distinct from Australian isolates, indicating local environmental, climatic or host features, rather than fungal features, govern the high incidence of A. flavus infection in this region. The results of this study have important implications for biological control strategies that aim to reduce aflatoxin by the introduction of non-toxigenic strains, as potentially any strain of A. flavus, and closely related species like A. minisclerotigenes, might be capable of human infection.

Draft Genome Sequence of Bacillus pumilus Fairview, an Isolate Recovered from a Microbial Methanogenic Enrichment of Coal Seam Gas Formation Water from Queensland, Australia.

Despite its global abundance, Bacillus pumilus is poorly studied. The Fairview strain was obtained from a methanogenic anaerobic coal digester. The draft genome sequence was 3.8 Mbp long and contained 3,890 protein-coding genes. Like the SAFR-032 strain, it includes B. pumilus-specific proteins that likely confer enhanced resistance to environmental stresses.

Draft Genome Sequence of Methanobacterium sp. Maddingley, Reconstructed from Metagenomic Sequencing of a Methanogenic Microbial Consortium Enriched from Coal-Seam Gas Formation Water.

The draft genome of Methanobacterium sp. Maddingley was reconstructed from metagenomic sequencing of a methanogenic microbial consortium enriched from coal-seam gas formation water. It is a hydrogenotrophic methanogen predicted to grow using hydrogen and carbon dioxide.

Draft Genome Sequence of Clostridium sp. Maddingley, Isolated from Coal-Seam Gas Formation Water.

Clostridium sp. Maddingley was isolated as an axenic culture from a brown coal-seam formation water sample collected from Victoria, Australia. It lacks the solventogenesis genes found in closely related clostridial strains. Metabolic reconstructions suggest that volatile fatty acids are the main fermentation end products.

Clostridium sporogenes PA 3679 and its uses in the derivation of thermal processing schedules for low-acid shelf-stable foods and as a research model for proteolytic Clostridium botulinum.

The putrefactive anaerobe Clostridium sporogenes PA 3679 has been widely used as a nontoxigenic surrogate for proteolytic Clostridium botulinum in the validation of thermal processes for low-acid shelf-stable foods, as a target organism in the derivation of thermal processes that reduce the risk of spoilage of such foods to an acceptable level, and as a research model for proteolytic strains of C. botulinum. Despite the importance of this organism, our knowledge of it has remained fragmented. In this article we draw together the literature associated with PA 3679 and discuss the identity of this organism, the phylogenetic relationships that exist between PA 3679 and various strains of C. sporogenes and proteolytic C. botulinum, the heat resistance characteristics of PA 3679, the advantages and limitations associated with its use in the derivation of thermal processing schedules, and the knowledge gaps and opportunities that exist with regard to its use as a research model for proteolytic C. botulinum. Phylogenetic analysis reviewed here suggests that PA 3679 is more closely related to various strains of proteolytic C. botulinum than to selected strains, including the type strain, of C. sporogenes. Even though PA 3679 is demonstrably nontoxigenic, the genetic basis of this nontoxigenic status remains to be elucidated, and the genetic sequence of this microorganism appears to be the key knowledge gap remaining to be filled. Our comprehensive review of comparative heat resistance data gathered for PA 3679 and proteolytic strains of C. botulinum over the past 100 years supports the practice of using PA 3679 as a (typically fail-safe) thermal processing surrogate for proteolytic C. botulinum.

Draft genome sequence of Clostridium sporogenes PA 3679, the common nontoxigenic surrogate for proteolytic Clostridium botulinum.

Clostridium sporogenes PA 3679 is widely used as a nontoxigenic surrogate for proteolytic strains of Clostridium botulinum in the derivation and validation of thermal processes in food. Here we report the draft assembly and annotation of the C. sporogenes PA 3679 genome. Preliminary analysis demonstrates a high degree of relatedness between C. sporogenes PA 3679 and sequenced strains of proteolytic C. botulinum.

Utility of microsatellite markers and amplified fragment length polymorphism in the study of potentially ochratoxigenic black aspergilli.

Microsatellite markers and the results of amplified fragment length polymorphism (AFLP) were compared in the characterization of 68 Aspergillus carbonarius and A. niger aggregate strains of differing ochratoxin-producing ability and from different geographic areas, isolated mainly from grapes and soil. AFLP was applied to both A. carbonarius and A. niger aggregate strains, and it clearly differentiated these species. Microsatellite markers were only applied to A. niger aggregate strains because of the species-specific nature of these markers. Both AFLP and microsatellite marker analyses were able to divide A. niger aggregate strains into the two recognized internal transcribed spacer (ITS)-5.8S rDNA RFLP types, N and T. Clustering of A. niger aggregate strains was similar in both AFLP and microsatellite analyses, yielding an additional separation of N type strains into two groups. Both microsatellite marker and AFLP analyses showed high levels of polymorphism in the A. niger aggregate (index of discriminatory power 0.991 and 1.0, respectively). Of the two techniques, microsatellite marker analysis was quicker and more straightforward to perform. In addition, microsatellite marker analysis is more reproducible, and the results can be expressed as quantitative data, making microsatellite markers a good candidate for use in large-scale studies of genetic diversity in A. niger aggregate species.

Xanthones from a microfungus of the genus Xylaria.

Chemical investigations of a microfungus Xylaria sp. isolated from the Australian rainforest tree Glochidion ferdinandi have afforded two new natural products, 2-hydroxy-6-methyl-8-methoxy-9-oxo-9H-xanthene-1-carboxylic acid (1) and 2-hydroxy-6-hydroxymethyl-8-methoxy-9-oxo-9H-xanthene-1-carboxylic acid (2). Compound 1 has previously been synthesised but only partially characterised. Methylation of 1 using diazomethane afforded the crystalline compound 2,8-dimethoxy-6-methyl-9-oxo-9H-xanthene-1-carboxylic acid methyl ester (3), whose structure was determined by single crystal X-ray analysis. This paper reports the full spectroscopic characterisation of compounds 1-3 by NMR, UV, IR and MS data. All compounds were inactive in a brine shrimp lethality assay and several antimicrobial screens.

Complete genomic sequence of SfV, a serotype-converting temperate bacteriophage of Shigella flexneri.

Bacteriophage SfV is a temperate serotype-converting phage of Shigella flexneri. SfV encodes the factors involved in type V O-antigen modification, and the serotype conversion and integration-excision modules of the phage have been isolated and characterized. We now report on the complete sequence of the SfV genome (37,074 bp). A total of 53 open reading frames were predicted from the nucleotide sequence, and analysis of the corresponding proteins was used to construct a functional map. The general organization of the genes in the SfV genome is similar to that of bacteriophage lambda, and numerous features of the sequence are described. The superinfection immunity system of SfV includes a lambda-like repression system and a P4-like transcription termination mechanism. Sequence analysis also suggests that SfV encodes multiple DNA methylases, and experiments confirmed that orf-41 encodes a Dam methylase. Studies conducted to determine if the phage-encoded methylase confers host DNA methylation showed that the two S. flexneri strains analyzed encode their own Dam methylase. Restriction mapping and sequence analysis revealed that the phage genome has cos sites at the termini. The tail assembly and structural genes of SfV show homology to those of phage Mu and Mu-like prophages in the genome of Escherichia coli O157:H7 and Haemophilus influenzae. Significant homology (30% of the genome in total) between sections of the early, regulatory, and structural regions of the SfV genome and the e14 and KpLE1 prophages in the E. coli K-12 genome were noted, suggesting that these three phages have common evolutionary origins.